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Abstract 

In supersymmetric theories, the Peccei-Quinn symmetry has a 
complex extension as a symmetry of the superpotential, so that the 
scalar potential always has an almost flat direction, the dilaton. We 
discuss how coherent oscillation of the dilaton affects axion cosmol- 
ogy. We stress that the dilaton decay, if its dominant mode is not into 
axions, releases large entropy at a late epoch of the Universe's evolu- 
tion to dilute axion energy density and the upperbound of the decay 
constant is raised up to about 10 15 GeV. The case of the M-theory 
axion is also discussed. 
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Although the standard model describes interactions of elementary parti- 
cles remarkably well, its extension is demanded by some fine tuning problems 
in it. One of such fine tuning problems is the strong CP problem QTJ| and 
another is the Higgs mass hierarchy problem. So far, the most attractive so- 
lution to the latter problem is supersymmetry 0, whereas the Peccei-Quinn 
(PQ) mechanism solves the former one in a beautiful way ||. The standard 
arguments based on astrophysics and cosmology constrain the PQ symmetry 
breaking scale to lie between 10 9 GeV and 10 13 GeV ||. |5|. 

In this paper, we shall consider the PQ mechanism in the framework of 
supersymmetry. Then the axion forms a supermultiplet. We shall discuss 
cosmological effects of the scalar superpartner of the axion and show that it 
can drastically change the standard axion cosmology. 

In supersymmetry, the PQ U{1)pq symmetry is extended to its complex 
form U(1) C P Q as a symmetry of the superpotential ||. To illustrate this point, 
let us consider a linear realization of the U(1)pq symmetry, in which a field 
transforms as 

<P* -> e iQi£ >\ (1) 

where a is a real parameter of the U(1)pq transformation and Qi the U(1)pq 
charge of the field. The point is that the superpotential W(4>), a holomor- 
phic function of the complex fields J , is invariant under the U(1)pq trans- 
formation with a being now an arbitrary complex number. When a is pure 
imaginary, the transformation corresponds to a dilatational one. Under the 
assumption of unbroken supersymmetry, this U(1) c P q symmetry leads to the 
existence of a non-compact flat direction in the scalar potential, associated 
with this dilatational transformation. We shall call, in this paper, the field 
along this flat direction a "dilaton" field, which will be denoted by (f>. Super- 
symmetry breaking generates the potential for the dilaton field and thus the 
dilaton mass. In the gravity-mediated supersymmetry breaking scenario, the 
dilaton mass, m^, is of the order of the gravitino mass, m 3 / 2 , which should 
be in the TeV range. 

Existence of the weakly-interacting scalar field, the dilaton, can drasti- 
cally change axion cosmology. In the early Universe the coherent mode of 
the dilaton field is in general displaced from its true minimum. As the ex- 
pansion rate of the Universe becomes comparable to the dilaton mass m^, 
the field starts damped oscillation. As we will see, the coherent oscillation of 
the dilaton and its subsequent decay may play important roles in the history 
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of the Universe. 

Cosmological implications of the dilaton were discussed in Ref. for 
a specific case where the PQ scale is radiatively generated. In the present 
paper, we stress that the existence of the flat direction is a general feature 
and explore its cosmological consequences in more general situations. 

We begin by discussing the decay rate of the dilaton. The dilaton gener- 
ally has a coupling to the axion through 

C eff = J-cfWadpa, (2) 

fa 

where F a is the U(1)pq symmetry breaking scale (or the axion decay con- 
stant) and / is a dimensionless parameter which depends on a model 
Generically it is of order unity but one can also construct a model with small 
/. From Eq. (0), it follows that the partial decay width is 

TU -> aa) = 1 1 (3) 

Besides the decay into two axions, the dilaton also decays to standard- 
model particles. First we consider a DFSZ axion model ||, in which the 
axion couples to quarks in the standard model at the tree level, whereas the 
coupling to gluons is suppressed by a one-loop factor. Since the coupling to 
a quark is proportional to its mass, the axion couples most strongly to the 
third-generation quarks. Explicitly it has the following coupling 

a [ 2x - 2x~ x - 1 

C eff = i— I — ■ T m t tj 5 t H ■ rm b bj 5 b + ...}, (4) 

F a [x + x~ v x + x^ 1 J 

where x = 1/ tan (3 is the ratio of the vacuum expectation values of the two 
Higgs scalar fields. The coupling of the dilaton is given by 

(p \ 2x - 2x~ 1 - 1 

C ef f = — { — ■ T m t tt H ■ -m b bb + . . . \ (5) 

F a [x + x~ L x + x _i J 

as a consequence of supersymmetry. Assuming that the channel to the it 
pair is open, one can calculate the decay rate of the dilaton 

3 f fm t \ 2 ( 2x \ 2 / 4m, 2 \ 3/2 
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Secondly in the case of a hadronic axion model |I(J, the coupling to the 
standard-model quarks vanishes at tree level. Among the standard model 
particles, the axion dominantly couples to the gluons at the one-loop order. 
The interaction is given by 

£e ff = g| (jr/^ + JrJ^J^) , (7) 

where g 3 is the SU (3)c gauge coupling constant and C is a model dependent 
constant of order unity. This yields the partial decay width into the two 
gluons 

We now consider cosmological effects of the decay of the dilaton. When 
the dimensionless parameter / in Eq. (0) is of order unity, the decay mode 
into the two axions tends to dominate over the other modes [§]. The standard 
big-bang nucleosynthesis (BBN) gives the upperbound of the energy density 
of the axions and the dilatons at the time of the BBN. Conservatively it must 
be less than the energy density contributed by one neutrino species. 

The dilaton field starts damped coherent oscillation when the expansion 
rate of the Universe becomes equal to the dilaton mass. The oscillating 
dilaton eventually decays into the axion pair. The produced axions do not 
get thermalized, but just red-shifted. When F a ^j 10 15 GeV, the life-time 
of the dilaton is shorter than 1 sec, and thus it decays before the BBN 
starts. If the energy density of the dilaton oscillation dominates the energy 
density of the Universe at its decay, the produced axions then dominate the 
energy density. This obviously upsets the standard BBN scenario and thus 
the dilaton-oscillation energy must not dominate the energy of the Universe 
before its decay. Therefore we will consider the case where the Universe is 
radiation dominated during the coherent oscillation. Then the ratio of the 
dilaton to the radiation energy density at the start of the oscillation is related 
to the same ratio at the time of the dilaton decay in the following way: 




b decay 



' osc 



(9) 



where t osc is the time when the oscillation starts and td eC ay is the dilaton 
life-time. If we write the initial amplitude of the dilaton field as eF a with e 
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being a dimensionless parameter, the above reads 



decay 



M 2 \ y V P 



(10) 



Here M represents the reduced Planck mass 2.4 x 10 18 GeV. The parame- 
ter e depends on the form of the scalar potential during inflationary epoch. 
Generically it is of order unity, but one can consider also the case where 
e ~ M/F a , i.e. the initial amplitude is of the order of the Planck scale. 
In the approximation of the instantaneous decay of the dilaton, the dilaton 
energy density is just converted to the produced axion energy density. We 
can evaluate the energy density of the produced axions p a at the time of the 
primordial nucleosynthesis (with the temperature ~ 1 MeV) as follows: 

(Pa) (pA ( g*(T D ) V /3 /32jT Fy ( gjT») \ 1/3 

UJlMeV \Pr) decay J V P M*™* \g.(MeV) J ' 

(11) 

where T D denotes the temperature at the time of dilaton decay. 

The observation of the primordial abundance of 4 He constrains the energy 
density of the produced axions at the BBN epoch to be less than that of the 
one neutrino species, which in turn means 

(12) 




The above constraint leads to a bound on the axion decay constant 

F a < 0.6 x 10 13 GeV ■ f^e^ (^) ^ (13) 

It is interesting to note that, for /, e =0(1), the upperbound of the axion 
decay constant fll~3|) coincides with the upperbound coming from the closure 
limit of the axion energy density due to misalignment. Note that in this 
region of F a the decay of the dilaton occurs before the nucleosynthesis and 
our analysis here is self-consistent. On the other hand, if e is much larger 
than unity, the above constraint is much severer than the conventional cos- 
mological upperbound. 
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When the axion decay constant exceeds 10 15 GeV, the dilaton still oscil- 
lates at 1 sec. At this moment, the ratio of the dilaton energy density to the 
radiation energy density is evaluated as 



The 4 He abundance constrains the above to be much smaller than unity. 
This in turn requires very small e, which is very unlikely. 

If the coupling of the dilaton to the axions, /, is small and the decay mode 
into the axion pair is not dominant, the upperbound 10 13 GeV for the PQ 
scale can be relaxed in an interesting way. Namely if the coherent oscillation 
of the dilaton field dominates the energy density of the Universe, then its 
decay followed by thermalization produces entropy, which may dilute the 
energy density of the axion oscillation. We shall look into this mechanism in 
detail. 

From the requirement that the decay of the dilaton does not spoil the BBN 
scenario, the reheat temperature after the entropy production, Tr, must be 
higher than 1 MeV. On the other hand, since the axion coherent oscillation 
starts around the QCD phase transition, the reheat temperature should be 
less than about 100 MeV in order that the entropy release can dilute the 
energy density of the axion oscillation. Because Tr is related to the decay 
width as 



the requirements on the reheat temperature constrain the axion decay con- 
stant, in which the dilution mechanism may work. 

We shall first consider the DFSZ axion. In this case, the dominant decay 
mode is into the top pair. Note that from Eqs. @ and (pj), we can evaluate 
the branching ratio as 



' 12 \m t J 

The BBN constraint requires that this should be smaller than 7/43. Thus / 
should be restricted to 




(14) 




(15) 




(16) 




(17) 
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From Eqs. (|6]) and (^) it follows that 

where we have neglected the decay mode to the bottom pair. From this, the 
requirement 1 MeV ^ 100 MeV leads to a constraint on the 
axion decay constant: 



m 1 ' 2 - « i9 » 

In order for the dilution mechanism to work efficiently, the energy density 
of the dilaton oscillation should dominate over that of the radiation already 
when the axion field starts oscillation. The ratio of the dilaton to radiation 
energy density at a temperature T is easily calculated to be 

E± „ (El) Z°££ (20) 

Pr \PrJ osc T 

where T osc is the temperature when the dilaton oscillation begins. From this 
we find that the temperature at which the dilaton energy equals the radiation 
energy is roughly 

r - 100MeV (lo^v) 2 (l^) 1/2 - < 21 > 

For the efficient dilution, this should be of the order of 100 MeV or more, 
which leads to 

eF^lO^GeV^)" 174 . (22) 

Keeping the above conditions in mind, let us next consider the cosmic 
abundances of the axion. As we discussed, if Tr^ 100 MeV, the reheating 
occurs after the axion field starts coherent oscillation and thus the entropy 
production dilutes the axion abundance. The abundance has been calculated 
in Ref . f 1] , with the result 

,2 _ . K ( T R \ ( F ad 



Q a h 2 - 5.3 * , (23) 

UMeVy \ 10 16 GeV/ ' V ; 
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where f2 a is the contribution of the axion energy density to the density pa- 
rameter, h is the Hubble constant in units of 100 km/s/Mpc, and 9 stands 
for the misalignment of the axion field (—71 < 6 < it). Plugging Eq. fllBp into 
the above equation, we obtain 

^ 2 - - 05 " 2 (lo4v)(l^) 1/2 - < 24 > 

From this equation we can conclude that the axion energy density does not 
overdose the Universe for 10 13 GeV ^ F a ^ 10 15 GeV, corresponding to the 
reheat temperature 100 MeV ^ Tr^ 1 MeV (see Eq. (0)), and thus the 
entropy production coming from the dilaton decay revives the DFSZ axion 
with the decay constant in the above rangeQ. Moreover, with F a ~ 10 15 GeV 
and the maximal value of 9 ~ tt, the axion will be able to constitute the 
dominant component of the dark matter of the Universe. 

On the other hand, if F a ^ 10 12 GeV and hence T R ^ 1 GeV, standard 
computation for the axion relic abundance should apply, namely 

F , 1.18 



n " h ~°- 23e (wtkv) ■ (25) 

This implies that the axion with F a ~ 10 12 GeV also provides a candidate 
for the dark matter. 

Next we turn to the case of the hadronic axion. Comparing Eq. (|3]) with 
Eq. (§), one finds that the dilaton dominantly decays into gluons, when 

f< ^ = O(10- 2 ) (26) 
for C = 0(1). When it is achieved, the reheat temperature is evaluated to 

With C = 0(1), the requirement Tr£ 1 MeV reads 

F a <10 14 Gev(^) 3/2 (28) 



lr The condition Eq. (||) is satisfied for e> O(l). 
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Note that when is around 1 TeV, the bound in the hadronic axion case 
is one order of magnitude lower than that in the DFSZ axion. 

For F a ~ 10 12 GeV-10 14 GeV, the reheat temperature is about 1-100 
MeV and therefore the entropy production by the dilaton decay dilutes the 
axion relic abundance^. Using Eqs. (|23|) and ([27D, one finds 

Q a h 2 ~ 5 x 10- 4 9 2 ( f a ) ( ^±-\ V2 . (29) 
Vl0 14 GeV/ VlTeV/ V ; 

This shows that the dilution mechanism makes the hadronic axion with F a ~ 
10 12 — 10 14 GeV cosmologically viable. Note, however, even with the maximal 
6 ~ 7r, Q a h 2 is at most 5 x 10 -3 , too small to be a dark matter candidate. 

When the U(1)pq breaking scale is as low as 10 11 GeV, the reheat temper- 
ature after the dilaton decay is 1 GeV. In this case, as was discussed for the 
other model, the dilaton decay does not dilute the axion energy density and 
hence Eq. ( p5|) can be used to estimate its relic abundance. For F a ~ 10 11 
GeV, the hadronic axion can marginally be the dark matter of the Universe 
if the maximum 9 ~ ir is taken. 

Finally we would like to discuss implications of our analysis to the M- 
theory axion. It has been argued ]I2|, 14] that the bulk moduli fields, liv- 



ing in the eleven-dimensional bulk, provide axion candidates in the massless 
spectrum of the strongly-coupled heterotic string theory (M-theory) com- 
pactified on a Calabi-Yau manifold. If the volume of the six-dimensional 
manifold is sufficiently large, the contributions of any kind of high energy 
origins to the axion potential will be strongly suppressed and the dominant 
contribution will come from the QCD anomaly. If this is the case the axion 
solves the strong CP problem. It can be shown that the coupling of the axion 
from the bulk moduli to other fields is suppressed by the (four-dimensional) 
Planck mass. This means that the decay constant of the axion is about 10 16 
GeV ||15|| . This axion is associated with its scalar superpartner whose mass 
is of the order of the gravitino mass. A study of the gaugino condensation 
scenario in the strong coupling case shows that the soft masses of gauginos 



and scalars are comparable to the gravitino mass [IB], and thus it should be 
in the weak scale to solve the naturalness problem on the Higgs mass. Then 
we immediately encounter a difficulty that the scalar partner of the axion 



2 For F a ~ 10 12 GeV, Eq. (^) requires a large initial amplitude of the dilaton field 
(e^ 0(10)) to achieve efficient dilution. 
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decays after the BBN commences (see Eqs. (0) and (pTD), with huge entropy 



produced, which is obviously a disaster |17| . To cure this, we have to invoke 
another entropy production mechanism with a higher reheat temperature 
(Tr^ 1 MeV) to dilute the scalar's energy density. One might imagine that 
faster decays of heavier moduli fields, if exist, could play this role. However, 
it is not the case, because these heavier fields start damped oscillation earlier 
than the dilaton does and the energy density is eventually dominated by the 
oscillation energy of the problematic scalar field (dilaton). To our knowledge, 



the only mechanism which can work is the thermal inflation ||18|| . This mech- 
anism can efficiently dilute the scalar's energy density to a harmless levelQ. 
Then the only constraint on the F a comes from the relic abundance of the 
axion. Applying Eq. to the present case, we conclude that, with the 
lowest allowed reheat temperature Tr ~ 1 MeV, the M-theory axion with 
F a ~ 10 16 GeV will survive the overclosure constraint if 9 is smaller than 0.3. 

In addition to the bulk moduli, there may exist moduli fields which live on 
the Eq boundary (boundary moduli). The boundary moduli has a coupling 
suppressed only by the GUT scale, a few times 10 16 GeV. If non-perturbative 
corrections of the superpotential of the boundary moduli are suppressed, it 
will be a dominant component of the QCD axion. Its scalar superpartner 
will get a mass comparable to the other scalars on the same boundary, e.g. 
sleptons and squarks. Applying again the argument of Ref. |15|] to estimate 



the axion decay constant, one infers F a ~ 10 13 — 10 14 GeV JT^|. Repeating 
our arguments for the hadronic axion case (see Eq. (|29"D ), we may conclude 
that the axion originated from the boundary moduli would be cosmologically 
viable by itself. However, the existence of the bulk moduli which seems 
to be inevitable again requires the late-time entropy production due to the 
thermal inflation. Then the QCD axion energy density will be totally diluted 
to a cosmologically negligible level. On the other hand, the other linear 
combination of the bulk axion and the boundary axion, which does not couple 
to QCD anomaly, will gain a mass from non-perturbative effects of string 
theory, such as a world-sheet instanton effect. It will be possible that it 
becomes the dark matter of the Universe after the dilution by the thermal 
inflation. 

To conclude, we have discussed the cosmological implications of the scalar 



3 Note that the reheat temperature should be lower than 100 MeV: otherwise the energy 
density of the axion oscillation would not be diluted. 
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superpartner of the axion, called the dilaton, in the gravity mediated super- 
symmetry breaking scenario. When the dilaton field dominantly decays into 
two axions, we have obtained the constraint on the axion decay constant 
from the BBN. On the other hand, when the dilaton dominantly decays into 
standard model particles, the decay can release entropy, which dilutes the 
energy density of the axion coherent oscillation. As a result the cosmologi- 
cal upperbound on the axion decay constraint put by the closure density of 
the Universe can be lifted up to about 10 15 GeV for the DFSZ axion and 
10 14 GeV for the hadronic axion, respectively. In the case of the M-theory 
axion, we have argued that the thermal inflation with low reheat temper- 
ature is necessary to dilute the energy density of the axion and its scalar 
superpartner. 
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